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Abstract 

We study the acceleration of the universe as a consequence of the time evolution of the vacuum energy in 
cosmological models based in braneworld theories in 5D. A variable vacuum energy may appear if the size 
of the extra dimension changes during the evolution of the universe. In this scenario the acceleration of the 
universe is related not only to the variation of the cosmological term, but also to the time evolution of G and, 
possibly, to the variation of other fundamental "constants" as well. This is because the expansion rate of the 
extra dimension appears in different contexts, notably in expressions concerning the variation of rest mass and 
electric charge. We concentrate our attention on spatially- fiat, homogeneous and isotropic, brane-universes where 
the matter density decreases as an inverse power of the scale factor, similar (but at different rate) to the power 
law in FRW-universes of general relativity. We show that these braneworld cosmologies are consistent with the 
observed accelerating universe and other observational requirements. In particular, G becomes constant and 
A(4) « const X asymptotically in time. Another important feature is that the models contain no "adjustable" 
parameters. All the quantities, even the five-dimensional ones, can be evaluated by means of measurements in 
4_D. We provide precise constrains on the cosmological parameters and demonstrate that the "effective" equation 
of state of the universe can, in principle, be determined by measurements of the deceleration parameter alone. 
We give an explicit expression relating the density parameters flp, Qa and the deceleration parameter q. These 
results constitute concrete predictions that may help in observations for an experimental/observational test of the 
model. 
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1 Introduction 



A number of recent observations suggest the possibility that our universe is speeding up its expansion pP-0]. The 
ratio of the (baryonic and cold non-baryonic dark) matter density to the critical density is approximately 1/3, while 
the remaining 2/3 comes from some kind of matter whose composition is unknown and which remains unclustered 
on all scales where gravitational clustering of ordinary (baryonic plus dark) matter is seen. Measurements of the 
luminosity- redshift relation observed for distant supernovae of type la, combined with data from the cosmic microwave 
background radiation, suggest that this new form of matter, which is sometimes called missing or dark energy, should 
possess a large negative pressure 

The simplest candidate for the missing energy is the cosmological constant [HI- Eli which looks like an ideal 
fluid with negative pressure pA = —p\. Another possibility is quintessence, which is a slowly-varying, spatially 
inhomogeneous field Ej-El whose effective equation of state pq = wqpq can vary with time and can take any value 
WQ > — 1. Alternative explanations for the acceleration of the universe, beyond dark energy, include phantom energy 
S2|-Qn]i certain modifications of general relativity ^El-EOli the gravitational leakage into extra dimensions |21)-|23j. 
Chaplygin gas |23-|2E1 as well as Cardasian models |2Z1-|2H1- 

In the present work, we study the acceleration of the universe as a consequence of the time evolution of the 
vacuum energy in cosmological models based on braneworld theories in 5D. Our aim is to examine the interplay 
between these theoretical models and observations, with some detail. Firstly, under very general conditions, we test 
whether these models are consistent with present cosmological observations. Secondly, we derive explicit expressions 
among the observational quantities and obtain the specific constraints that these models impose on the cosmological 
parameters. 

In brane theories our four-dimensional universe is depicted as a singular hypersurface |29| . a 3-brane, embedded 
in a 4 -|- d dimensional spacetime, or bulk. The scenario is that the matter fields are confined to the brane, while 
gravity propagates in the extra d dimensions as well. The effective equations for gravity in AD were obtained by 
Shiromizu et al |3()j by using Israel's boundary conditions and imposing Z2 symmetry for the bulk spacetime, about 
our brane- universe fixed at some y = yo- These equations predict five-dimensional corrections {d — 1) to the usual 
general relativity in AD and, what is specially relevant to us, they provide explicit expressions that link the tension 
of the brane to the observed gravitational coupling G and cosmological term A(4) . 

Most of the studies of braneworld cosmologies in 5D deal with the case where the size of the fifth dimension 
is constant, i.e., 344 = const This is consistent with the assumption of time-independent vacuum energy, 

which in turn leads to that A(4) and G be constants automatically. Thus, if the dark energy is identified with a 
time-independent vacuum energy (positive cosmological constant), then the universe will expand forever, the energy 
density of all kinds of matter, except for vacuum, tend to zero exponentially, and the spacetime metric approaches 
the de Sitter one. 

However, this simple picture is altered if the vacuum energy is allowed to vary with time. The crucial point here is 
that the introduction of a time- varying 544 in braneworld theory allows the construction of a number of cosmological 
models in AD which have good physical properties but do not admit constant vacuum energy |33| . This physical 
scenario cannot be ruled out a priori. After all, the constancy of 1744 is an external condition and not a requirement 
of the field equations. In this scenario A(4) and G vary in cosmological time, which in turn influences the expansion 
of the universe, through the generalized Friedmann equation. 

In a recent work |34| we studied the cosmological consequences of a braneworld theory, with time dependent 
vacuum energy, where the time variation of G can be expressed as {G/G) — gH with \g\ < 0.1, in order to satisfy 
observational requirements. The main results emerging from that work were that: (1) Our universe must be spatially- 
flat {k = 0), (2) the models with zero and negative bulk cosmological constants agree with the observed accelerating 
universe, while fitting simultaneously the observational data for the density and deceleration parameter, and (3) 
the universe must recoUapse if it is embedded in an anti-de Sitter bulk, although we are nowhere near the time of 
recoUapse. 

The question arises of whether or not these results are typical in braneworld cosmologies with variable vacuum 
energy. In order to elucidate this question we need to examine other models of this kind, based on the simplest and 
most general assumptions possible. In this work we use two assumptions that are common in modern cosmology: 

1. The matter density decreases as an inverse power of the scale factor. 
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2. The universe is spatially flat. 

Our aim in this work is to study the cosmological consequences of the above assumptions on the evolution of 
braneworld cosmologies. We will discuss the following questions: What is the general evolution of the scale factor? 
Does it agree with the observed accelerating universe? Is the variation of G and A(4') compatible with observations? 
What is the effective equation of state of the universe? What is the most general relationship between density and 
deceleration parameters? Based on our model and observational data, how to estimate and/or constrain the value 
of cosmological parameters?. Although the model is classical, we also tackle the question of whether the vacuum 
energy in our model can be identified with the energy of some scalar field rolling-down its potential. 

The main conclusions from our work are the following: 

The general behavior of the braneworld cosmologies under consideration can be separated into three stages. The 
first one is the early universe, when the quadratic correction term is dominant and the deceleration parameter is 
positive (q > 0). As the universe expands and ages, the quadratic term decreases while the terms for matter and 
vacuum increase. The second stage in when the vacuum energy density and the matter energy density become equal 
to each other. This occurs for fip — 1/2. At this stage the effects due to vacuum become important; the matter 
density parameter flp starts decreasing (instead of increasing) and q approaches zero. In the third stage the universe 
is dominated by the vacuum energy, which is responsible for the observed present acceleration. Since the quadratic 
term is now negligible, the universe can be considered as a mixture of dust and a variable cosmological "constant" . 
However, the effective equation of state of the universe, for the total pressure and density, is quintessence-like namely, 
pt Ri — 0.7pf, which coincides with the one obtained from other models based on FRW cosmologies |r{5 | -|36 | . 

From a conceptual viewpoint, these cosmologies present a consistent picture where the current acceleration of the 
universe is not an isolated phenomenon, but it is intrinsically interwoven with the time variation of the cosmological 
parameters G and A(4). Since the expansion rate of the extra dimension appears in different contexts such as the 
variation of rest mass |37| and the fine-structure-constant 38 , the braneworld scenario may provide us a theoretical 
framework to unify all these, apparently, separated phenomena as different consequences in AD of the time evolution 
of the extra dimension. This is a new step toward understanding how the universe works. 

From an observational viewpoint, based on the positivity of the mass/energy density and the observational fact 
that l^l < 0.1, we provide precise constrains on the cosmological parameters. We show that the equation of state of 
the universe can, in principle, be determined by measurements of the deceleration parameter alone. We obtain an 
explicit expression between the density parameters f2p, f2A and the deceleration parameter q. These results constitute 
concrete predictions that may help in observations for an experimental/observational test of the model. 

This paper is organized as follows. In section 2 we give a brief summary of the equations for homogeneous 
cosmologies in 5D based on braneworld theory. In section 3 we discuss the mathematical consequences of the 
two assumptions mentioned above. The integration of the equations is presented in section 4. Some of the physical 
consequences are explored in section 5; we derive relationships among the observational parameters, including the age 
of the universe and the effective equation of state of the universe. We also show how to evaluate the five-dimensional 
quantities that appear in the four-dimensional equations. In sections 6 we study the experimental/observational 
restrictions on the cosmological parameters. In section 7 we present an outline of the evolution of the universe using 
the results of our model. Section 8 is a summary. Finally, in the appendix we construct the scalar field potential 
corresponding to the vacuum energy in our model. 



2 Homogeneous cosmology in 5D 

In order to facilitate the discussion, and set the notation, wc start with a brief summary of the pertinent ideas and 
equations in the braneworld scenario. In this scenario our homogeneous and isotropic universe is envisioned as a 
singular hypersurface embedded in a five-dimensional manifold with metric 



(1 - fcr2) 



+ r^ide^+sm'^ 0d(l)^ 



^\t,y)dy\ 



(1) 
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where t, r, 6 and 4> are the usual coordinates for a spacetime with sphericaUy symmetric spatial sections and k = 
0, +1, — 1. The metric is a solution of the five-dimensional Einstein equations 

(5)G^B = ^"^Rab - IgAB^-'^R = klf^TAB. (2) 

where '^^^Tab is the five-dimensional energy- momentum tensor and /c(5) is a constant introduced for dimensional 
considerations. 

The energy-momentum tensor on the brane t^^, is separated in two parts, 



^9 1 



(3) 



where a is the tension of the branc in 51?, which is interpreted as the vacuum energy of the braneworld, and T^j^ 
represents the energy-momentum tensor of ordinary matter in AD. 

There are two assumptions relating the physics in 4D to the geometry of the bulk. The first one is that bulk space- 
time possesses Z2 symmetry about the brane. As a consequence the matter in AD becomes completely determined 
by the geometry in 51? through Israel's boundary conditions. Namely, for a perfect fiuid 



with 

the matter density is given by 
and 



P = 7P, 



(-2) 

bra7ie 



a n 
a n 



brane 



^(5)(7+l)$|br 



(37- 



a n 



bra 



where a prime denotes derivative with respect to the extra variable y. 

The second assumption is that the brane is embedded in an Anti-de Sitter bulk, i.e., ^^^T^^ is taken as 



(5) 



Tab = ^{5)9AB, 



(4) 
(5) 
(6) 

(7) 
(8) 



where A(5) < 0. Since ^^•'r^4 = 0, it follows that the energy-momentum tensor on the brane r^^ is a conserved 
quantity, viz., 

V% = 0. (9) 
Another important consequence of ||SJ) is that the field equations Q admit a first integral, namely, 



fc(5)A(5) 

6 



a 

a<i> 



C 



(10) 



where C is a constant of integration which arises from the projection of the Weyl curvature tensor of the bulk on 
the brane. Evaluating (|10|l at the brane, which is fixed at some y — ytrane = const, as well as using Q and Q, we 
obtain the generalized Friedmann equation, viz.. 



3(^ 

ao 



. o ^ (5) ? 3fc 3C 

_L ^ tig Q 



where ao(t) = a{t,ybrane), and 



1 / kf^^c 



(11) 



(12) 
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SttG^^. (13) 

These quantities are interpreted as the net cosmological term and gravitational couphng in 4 dimensions, respectively. 

Equation Hll() contains two novel features; it relates the fundamental quantities A(4) and G to the vacuum energy, 
and carries higher-dimensional modifications to the Friedmann-Robertson- Walker (FRW) cosmological models of 
general relativity. Namely, local quadratic corrections via , and the nonlocal corrections from the free gravitational 
field in the bulk, transmitted by the dark radiation term C/a*. 

Except for the condition that n = 1 at the brane, the generalized Friedmann equation is valid for arbitrary 
<i>(t, y) and n(i, y) in the bulk |32|. This equation allows us to examine the evolution of the brane without using any 
particular solution of the five-dimensional field equations. In what follows we will omit the subscript 0. 



3 Variable vacuum energy 

In equation Hll() . G and A(4) are usually assumed to be "truly" constants. However, the vacuum energy density a 
does not have to be a constant. From ((TJ it follows that it depends on the details of the model. Indeed, we have 
recently shown j33| that there are several models, with reasonable physical properties, for which a variable <& induces 
a variation in the vacuum energy a. 

For variable vacuum energy, the conservation equations for a perfect fluid which satisfies (|SJl, yield 

/9 + 3p(7-|- 1)- = -CT, (14) 
a 

For the case of constant a, we recover the familiar relationship between the matter energy density and the expansion 
factor a, viz., 

For a given ct as a function of a, a = cr(a), we integrate (|14|) and substitute the resulting function p = p{a) into l|ll|l . 
and thus obtain the corresponding Friedmann equation. 

The variation of the vacuum energy is deeply rooted in fundamental physics. The simplest microphysical model 
for a variable A(4), as well as for quintessence, is the energy associated with a slowly evolving cosmic scalar field 
<j) with some self-interaction potential V{(j)) minimally coupled to gravity. The potentials are suggested by particle 
physics, but in principle V{4>) can be determined from supernova observations |39|-|41j. 

In this work, instead of constructing a field theory model for the time evolution of the vacuum energy, we adopt 
a phenomenological approach. Namely, we assume that the energy density p can be expressed in way similar to l(TK|l . 
i.e., that it decays as a power of the scale factor a. Specifically, 

^ /5^-l, (16) 



a3(/3+i) 



where Z) is a positive constant and 3(/? -I- 1) is just a useful way of writing the power. This is convenient in order 
to simplify the comparison with the case where a = const, and with some previous results. With this assumption, 
equation 1)14(1 can be easily integrated as 

^ = ^"+(/3+l)a3(/3+i)' 

where ctq is a constant of integration. Thus, any deviation from the standard power law 1(15(1 . i.e., /3 7^ 7, yields a 
dynamical behavior for ct, implying that G and A(4) are functions of time. At this stage, we consider it premature to 
justify this model on the basis of fundamental physics, however in the appendix we show that 1(17(1 can be obtained 
from a simple potential V{(j)). Our main purpose now is to explore whether 1(17(1 leads to an scenario where the 
evolution of the scale factor, as well as the simultaneous variation of G and A(4) are consistent with cosmological 
observations. 
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With this aim we substitute Hlt)|) and (|17|l into Hll|l to obtain 



, a 
3 ( - 

a 



^(5) 

2 

U4 



A 



(5) 



^(5)^0 /7 + 1 



6 



D 



'^(5) 



7+1 



3fc 

Z2 



3C 



(18) 



This equation for /? = 7 reduces to the generalized Friedmann equation for constant vacuum energy. However, for 
/3 7^ 7 it generates models with varying G and A(4-) . 

In order to discuss such models, it is convenient to introduce the quantities 



in terms of which (|18|) becomes 



.3(/3+l) 



A 

B 
C 



k^ 

''(5) 



(7- 



(7+l)(/3 + l)i?, 



^(/5 + lf (A 




(19) 



l^j.{ef>+i)/{3fJ+3) _ ^^(6/3+2)7(3/3+3) 



(20) 



This equation admits exact integration, in terms of elementary functions, for a wide variety of parameters A, B, 
C, P and C. In this work we assume k = C = 0, which assures that at late times a ^ p our brane-model becomes 
indistinguishable from the flat FRW universes of general relativity. 

At this point let us notice that, the astrophysical data from BOOMERANG and WMAP indicate that our universe 
(which is assumed to be described by general relativity) is flat {k = 0) and A-dominated. Besides, the constant C, 
which is related to the bulk Weyl tensor and corresponds to an effective radiation term, is constrained to be small 
enough at the time of nucleosynthesis and it should be negligible today. Thus, setting fc = C = in (|20|) we assure 
that at late times our model is compatible with observations. 



4 Spatially flat universe 

In this section we discuss the solutions of (|2Uf) with k = C = 0. The time evolution of the brane-universe under 
consideration crucially depends on C. Namely, x ~ sinh ^/Ct for C > 0, and x ~ sin \J —Ct for C < 0. Let us notice 
that C can be written as C = 3(/3 + l)^Ae// with 



Ae// = '^k- A(,) + ^ . (21) 




This quantity may take any arbitrary value by appropriately specifying the constants A(5) and/or do- It plays the 
role of an "effective" cosmological constant^ in (|18|) . If it is positive the scale factor can take arbitrary large values. 
Eventually, it dominates the evolution causing an exponential expansion as in the de Sitter model. If it is negative, 
the scale factor is bounded above and ultimately the universe recoUapses. 

^In the case of constant vacuum energy 1211 is indeed the cosmological constant in 4Z). 
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In cosmological models with constant vacuum energy, the standard assumption is the vanishing of the net cosmo- 
logical term A(4) defined as in H12|) . This suppresses the exponential expansion and the recoUapse mentioned above. 
In the case of varying a the equivalent assumption is 



Ae// = 0. 

In order to avoid misunderstanding we should note that H22|l does not imply Ku\ — 0. It requires either A 



(22) 



fc(5)'^o/6 - or A(5) = -k' cj'jQ < 



H5) 



4.1 Bulk with A 



(5) 







For the first case, with (Tq = in 117|) . the scale factor is given by 



a 



"^^+'^ ^^h + l)D{t-to), A(5)=0, ao = 0. 



(23) 



where to is a constant of integration. The physical interpretation of this solution has already been discussed in our 
previous work 34 . It represents a brane-universe where G and A(4) vary as 



G 
G 



H, A 



(4) 



(24) 



and the ratio (p/o') = (/? + l)/(7 — P) remains constant in time, although varies with the epoch, i.e., with 7. In the 
remainder of this work, we will consider in detail the properties of the second case (A(5) = —k'^^^a^/G < 0). 

4.2 Anti de Sitter bulk 

In this case the expansion of the universe is governed by 



,3(/3+l 



)_ ^(7 + 1) 



fc'5)'^o(/?+l) (^-^o) -4 , A(5)<0, ao^O, 



Sao iP + 1) 

where to is a constant of integration. With this expression we find the "deceleration" parameter q = —da/a? as 

1 + 3/3 6 



q = 



iP+l)kf.,a^oit~tor' 



the Hubble parameter H = a/a 



the vacuum energy 



and the ratio 



1 

H 



3(/?+l) 



{t - to) 



a 

— = 1 



(/3 + 1)4)^2 (t-io)' 
8(7 - /3) 



(7 + 1) k\^^al{(i+\f(t-tof-^ 



In order to assure a(t) > and H > 0, for all t > 0, we choose the constant of integration to as 

2 



to 



k^^CToiP+l)- 



(25) 
(26) 

(27) 
(28) 

(29) 
(30) 



Besides, for this choice the big bang occurs at t = 0. Now the positivity of G from (|13|) demands cr > 0. On the 
other hand, selecting 7 > /?, we assure (ct/ctq) > for all values of t. Consequently ctq must be positive. 
In the remaining of this work we will study the physical properties of this solution. 
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5 Cosmological parameters 



In order to compare our model (|25|l - (|30|l with observational data we need to derive relations among the observational 
quantities. Although the generalized Friedmann equation 1)11(1 does not depend on the five-dimensional model, the 
task here is more complicated than in the FRW models because, we have to deal with the quantities fc(5), a. 

It is not diSicult to show, using the above formulae, that the density and deceleration parameters, Hp and g, are 
given by 

SttGp _ ^ _ 2(/3+l)[(/j + l) + (7-/3)(p/ao)] 

3i/2 P (7 + l)[2(/3+l) + (7 + l)(p/c7o)]' "^"^^ 

and 

^ 1 + 3/? 3(7 + l)(/3+l)(p/ao) . . 

2 ^2[2(/3+1) + (7 + 1)(pM)]- ^ ^ 

Thus, for any value of 7 (say dust 7 = 0), if we are able to measure ilp and g, then we obtain the present values of 
(p/cTo) and (5. These can be used to evaluate the cosmological "constant" A(4-) in terms of the Hubble parameter. 
Namely, 

. _ 3(7-/?)[2(/3+l) + (7-/3)(p/^o)] ,,2 zoos 
(7 + l)[2(/3+l) + (7 + l)(p/ao)] " ^''^ 

We note that A(4) here is not strictly proportional to iJ^ because (p/ctq) is a function of time. They become 
proportional only asymptotically, for (p/ctq) ^ 0. A similar situation happens with the ratio GjG where we have 

§ = -3(7-/3) (9 i^. (34) 

So, unlike in our previous study [3^, in the present work GjG ^ const x H. It is important to note that the 
five-dimensional quantities fc(5) and A(5) can be evaluated by means of quantities measured in 4D only. Specifically, 
for fc(5) we have 



4 _ ASttG 
«(5) - 



(7 - /3) / P 



(/3 + 1) V^o 



(35) 



P 

Similarly, from H21|l and (|22ll we obtain Aj-s). 

Let us now construct a working expression for the age of the universe T. First we isolate (t — to) from ()16|l and 

as 

_ , ^2 _ 4[2(/3 + l) + (7-H)(p/ao)] 
^ °^ -(p/ao)(7+l)(/3 + l)^(fcf5)'-§)' ^ ^ 



Next, we substitute this into from H27|l and obtain 

^(5)^0 _ 6(/3 + l) 



H 



y/{p/ao){j + miP+l) + (7 + l)(p/^o)] 



(37) 



This is a dimensionless quantity, which is perfectly well defined for each set {Hp, q). Finally, from H3U|) . H3()|l and (|37|l 
we obtain the desired expression, viz.. 



^ [2(/3 + l) + (7 + l)(p/ao)] [ ^ _ V(7 + l)(p/ao) 1 

3(/3 + l)' 1 v/[2(/3 + l) + (7 + l)(p^o)] J 

5.1 Parameters as functions of q and fip 

It is useful to express the above quantities solely in terms of the observables q and flp. To this end we use ((31|l and 
(132(1 to eliminate /3 and p/aa. 
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From H32|) we obtain p/cfq as a function of q and (3. Namely, 

p (/3 + l)[2(q + l)-3(/3+l)] 



(39) 



ao (7 + l)[3(/?+l)-(g + l)] ■ 

On the other hand, from (|31|l we obtain (p/ctq) as a function of [3 and fJp. From these two equations we get a 
quadratic equation for (/3 + 1). Namely^, 



6(/3 + 1)2 - A{q + l)(/3 + 1) + [2(7 + l){q + 1) - 3(7 + if^p] = 0. 
In order to solve this equation, is important to notice that the condition < (p/tro) < cxd in 1)39(1 requires 



1 



< 



13+1 
9+1 



< 



Consequently, the physical solution for /3, i.e. the one which leads to a positive p/tro, is given by 

g + 1 



(/? + i) 



3 



1 



'1-3 
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+ 2 







(40) 



(41) 



(42) 



The second solution, with a negative sign in front of the root, is outside of the allowed region (|41|) and leads to a 
negative p/cjq. 

In order to simplify the notation, let us introduce the auxiliary function 



/(g,n„7) = i + i 1-3 



2 Vg + 1 



and denote 



With this notation, we have 



7+1 ^ , 9/^7+1 
9 + 1 

g + 1 

7+1' 



(43) 
(44) 



(/3 + 1) 

P_ 

a 

TH 



(7 + 1), 



e/(2-/) 



3-e/(2-/)' 

^(3-e/)[3-e(2-/)] 



1 



P ^ ^/(2-/) 
fTo 3(/-l)' 



/2-_/ 



(45) 



(7 + l)(/-l)e 

In the above equations flA — A(4)/3if2 and 1 < f{q,V,p,j) < 2. The positivity of G requires (7 — /3) = (9 + 1)(3 — 
C/)/3^ > 0, which in turn imposes an upper bound on namely, 

e < 3. (46) 

It is interesting to note that G decreases in time. Namely, from (|34|) and the above expressions, we obtain 

G ^ e./(3-e./)(2-./)(l + 7) 
G 3-e/(2-/) 



-H. 



(47) 



Clearly, G/G < 0, as long as (7 + 1) > 0, which means that gravity was more important before than at the present 
time. 

In addition we note that any combination of q and 7, giving the same ^ in (|44|l . yields identical values for p/a 
and r^A- Only f3, the age of the universe T and G/G are sensible to the specific choice of 7 in the equation of state 
for ordinary matter. 

^Notice that Qp and q separate for (7 + 1) = 2(/3 + 1), i.e., /3 = (7 — l)/2. For this value, the density parameter is constant, viz.. 
Up = 1/2 for all q. 
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5.2 The equation of state of the universe 

In order to complete our set of equations, we find the equation of state for the total energy density pt and pressure 
Pt of the universe. It can be written as 

Ei^w, (48) 
Pt 

where 

^-^+[2(/? + l) + (7 + l)(p/ao)] - 3 ■ ^^'^ 
The value of w is not constant and differs from the equation of state ^ = p/p for ordinary matter on the brane. 
Thus, we arrive at an important result: the equation of state of the universe can be determined by measurements of 
the deceleration parameter alone. 

Finally, we would like to note the following relationship between the observational quantities flp, q and 7. 
Namely, 

2£i±i£[o,-i^^]+,(l-7)-7(l-3«-2^0^ (50) 

Since the overwhelming time of the evolution of the universe is spent in the matter-dominated {p — 0) domain, this 
equation is of practical use for 7 = . Thus, 

q^2-^inp + 2QA). (51) 

Consequently, for the equation of state of the universe we get 

w = l-np- 2rjA. (52) 

It is important to mention that the equations in this section reduce to the appropriate FRW-counterparts for J^a = 
and rip — 1. In particular we obtain q — 1/2 and w = as in the dust-FRW cosmologies. In the next section we will 
use the above expressions to evaluate the parameters of the model. 



6 Physical restrictions on {Qp, q) 



Thus, a wealth of important information about our universe can be derived from (^p,q). However, not all sets of 
values (ilp,?) can generate adequate physical models. The object of this section is three-fold. Firstly, to analyze 
the constraints that observation can place on them. Secondly, to show that our model agrees with the observed 
accelerating universe, and fits observational data. Thirdly, to estimate the accuracy with which Hp and q can be 
constrained by observations, based on our model. 

6.1 The positive-energy condition 

We have already used the condition p > to single out the physical solution of /3. But this is not all, the positive- 
energy condition is a powerful tool which allows us to discriminate between what is possible or not. Here we use this 
condition to formulate some specific constraints on q and/or flp. 

To this end, we substitute our auxiliary function f{q,rip,j) into (|in|) and obtain, 

ln.-n-.(i±l^)-.{i±i)'a,. (53, 

Then, from the positive-energy condition (|41|l . which in terms of / becomes 1 < / < 2, it follows that 

2 .o^7 + l\' 0^7+1 



This is an important expression because it involves observational quantities only, viz., q, Qp and 7. We now proceed 
to discuss its physical implications. 
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6.1.1 Restrictions on C = + l)/(7 + 1) 



We have already mentioned that for /? = (7 — l)/2 the density parameter is constant and equal to 1/2, regardless of 
the specific value of 7 and q. We can easily show that the solutions that become vacuum dominated are those with 
/? < (7 - l)/2. Indeed, from iglj and (EU we get 

f^A«^, (55) 
7+1 7+1 

asymptotically, as (p/co) ^ 0. In this limit, fl\ = flp for /3 = (7 — l)/2; the vacuum energy ultimately dominates 
only if /3 < (7 — l)/2. For these solutions flp is always positive and less that 1/2. 

The left-hand-side of (|54|l provides a quadratic expression for i^, which for fl < 1/2 requires either ^ < (3/2) [1 — 
^yi — 2^lp] or f > (3/2) [1 + \/l — 2rip], while the right-hand-side requires ^ > 3r2p/2. Thus, for every given value of 
51p < 1/2, there are two regions where the solution is consistent with p > 0. These are 

9 + 1 



and 



7 + 1/ 2 



31)p (q^l 



>-[i + v/r^2f7;], (56) 



2 <lTTTj<2[^-v/^^]- (57) 
We now proceed to discuss the physical properties of the solution in each region. 

Decelerated expansion: Let us calculate the cosmological parameters in the region controlled by (|56|l . If we 
take fip — 0.1, then the allowed values of ^ are spread over a small range, viz., 2.842 < ^ < 2.850. The lower limit 
comes from (|56f) . while the upper one comes from G > 0, in agreement with (|46|l . For these values, J^a is about 
0.001, indicating that the evolution here is dominated by the quadratic correction term f2p2. The huge gravitational 
attraction produced by this term (rip2 « 0.899) explains the large deceleration parameter q k, 1.846 + 2.8467, which 
is positive for all forms of matter including dust (7 = 0), radiation (7 = 1/3), stiff matter (7 — 1) and a network of 
cosmic strings (7 = — 1/3). 

The same picture emerges if we use other values for fip in H56|l . Consequently, if the observational parameters 
(rip, g, 7) are tied up as in (|56|l . then the solution is relevant to the description of an universe in decelerated expansion. 
We should mention at this point that there is plenty of observational evidence for a decelerated universe in the recent 
past, see e.g. 1121-113 ■ 



Accelerated expansion: A similar calculation, but in the regime governed by l|57(l . yields a totally different 
picture. Namely, the universe is dominated by the A(4)-term and is expanding with acceleration. In view of the 
relevance of this type of solutions to the recent discovery of an accelerating universe, we will examine the cosmological 
parameters in more detail. Our aim is to see whether they fit observational data. 

In table 1 we consider the behavior of dust-filled universes for several values'^ of fip and study how H57|l restraints 
the cosmological parameters. We use the mean value of q to calculate /? and then obtain the "present-day" value 
of {p/a)^ Ha, the age of the universe T, as well as other relevant cosmological parameters, such as w — Pt/pt and 
g = G/HG. 



Table 1: Determining the cosmological parameters from flp 



np 


q 


/3 




r!A 


HT 


w 


l5l 


0.3 


-0.499 + 0.051 


-0.729 


0.113 


0.683 


1.659 


-0.666 


0.248 


0.2 


-0.681 + 0.019 


-0.823 


0.063 


0.794 


2.600 


-0.787 


0.156 


0.1 


-0.846 + 0.004 


-0.913 


0.027 


0.899 


5.379 


-0.897 


0.073 



reliable and definitive determination of f2p has thus far eluded cosmologists. Currently, the observational data based on dynamical 
techniques indicate that Qp fx 0.1 — 0.3 seem to be the most probably options. 
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Table 1 shows that the universe is expanding with acceleration, in agreement with modern observations. The 
acceleration is driven by the repulsive effect of the "dark energy" associated with A(4-) , which clearly dominates the 
evolution here. For flp w 0.3 the quadratic correction is "negligible", flp2 « 0.02, and decreases as the universe gets 
older. 

Thus, at late times Qp + ^Ia « 1. In other words, the resulting cosmological model is a combination of ordinary 
matter with an unknown form of matter which is effectively described by a cosmological term which varies with time. 
The latter satisfies the equation of state pA — — Pa, with positive energy density, and dominates over the regular 
matter. As a consequence, at the present stage, the evolution of the universe is driven by an effective equation of 
state which is "quintessence-like", in the sense that (—1 < w < 0). 

We note that, knowing the value of the matter density p, the vacuum energy density can be obtained from the 
fourth column in table 1. We can also calculate the constant cto, or p/aQ, from l|45(l . This allows us to obtain A;(5) 
from H37|l as well as the value of the cosmological constant in 5D, A(5) = — /c^g-jCr^/G, which follows from H21|l and 



6.1.2 Restriction on flo 



Reciprocally, H54|l can be used to restrict the density parameter by measurements of q. Namely, for every given value 
of q, there is a range of possible values for ftp, viz.. 



2(q 
3 



7 + 1 



1 - 



111 

7 + 1 



<flp< 



g+1 
7 + 1 



(58) 



If we take the deceleration parameter as in the usual dust FRW cosmologies, i.e., q = +0.5, then 1/2 < i7 < 1. 
However, recent measurements indicate that, the current universe is speeding up its expansion with an acceleration 
parameter which is roughly g = — 0.5±0.2. If we fix q today, then the adequate values of flp are spread over a small 
range, which decreases as the universe ages. In table 2 we use the mean value of Hp to calculate the cosmological 
parameters, corresponding to various applicable values of q, for universes in the dust-dominated era (7 = 0). The 
density parameter flp can be determined from l|58() with an accuracy which is up to 86% — 95%. 



Table 2: Determining the cosmological parameters from q 



q 


rip 








HT 


w 


l5l 


-0.3 


0.412 + 0.054 


-0.602 


0.133 


0.561 


1.183 


-0.533 


0.240 


-0.5 


0.305 + 0.028 


-0.717 


0.093 


0.681 


1.657 


-0.666 


0.200 


-0.7 


0.190 + 0.010 


-0.829 


0.052 


0.805 


2.761 


-0.800 


0.131 



It is important to note that the possibly values of Sip {Vl w 0.1 — 0.3) and q {q — —0.5 + 0.2) contain an error, 
or uncertainty, of about 50% and 40%, respectively. This is a clear evidence of the difficulties involved in the 
determination of Qp and q from observations. We can use our model to narrow down these uncertainties and get Qp 
and q with a better accuracy. 

First, from table 1 we see that for SI w 0.1 — 0.3 the acceleration parameter varies in the range —0.85 < q < —0.448. 
Then, from the last row in table 2 we see that large negative values for q, bellow —0.7, are possible because we are 
considering density parameters less than S7p = 0.18. Based upon the observational limit on g, we can set a lower 
bound to Qp. Consequently, 



0.18-0.3. 



(59) 



Second, from table 2 we see that that q = —0.5 + 0.2 requires Sip ~ 0.18 — 0.446. Then, from the first row of table 1 
we see that large values for Sip, above 0.3, result from considering q greater than w —0.448. Again, based upon the 



observational limits on the density parameter, we can set q — —0.448 as the upper limit of q. Thus, 

« = -0.574+ 0.126. 



(60) 



The conclusion is that our model allows to significantly reduce the uncertainty from 50% to 25% for Sip, and from 
40% to 22% for q. We emphasize, however. Sip and q are not completely independent: fixing f2p today also fixes q 
with a 90 — 99% confidence. 
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6.2 Constraints from the variation of G 

We now proceed to show that our model is compatible with the observational constraints on the variation of G. To 
this end, we use H47|l and denote 

e/(3-e/)(2 -/)(! + 7) ...s 
^ 3-e/(2-/) ' ^^'^ 

so that G/G = gH. In astronomical and geophysical experiments, measurements of possible variations of G with 
cosmic time provide experimental bounds on G/G, which can be summarized as** 

l5l < 0.1. (62) 

Clearly, the positivity of energy is not enough to guarantee the fulfillment of this condition. This is illustrated in 
tables 1 and 2. Thus, H62|l imposes observational constraints on {Up,q). 

For example, if we take Op = 0.3 and 7 = 0, then H()2|) requires —0.571 < q < —0.529. On the other hand, 
(|57|l requires —0.55 < q < —0.448. As a consequence, the physical range of q reduces to —0.55 < q < —0.529, or 
q = —0.539 ±0.011. Therefore, for a fixed flp, the best lower bound on q is provided by the positive-energy condition 
(jsTjl . while the best upper bound is provided by (|^ . 

Table 3 shows how the observational evidence \g\ < 0.1 allows to increase the precision in the determination of 
q to be more than 98%. Again, we use the mean value of q to calculate the rest of the parameters. We note that 
smaller Op means that the universe is older and /? is closer to —1 today. We also note that 

w « — Oa, (63) 

with an error of less than 1%, which is consistent with H52|) for an "old" (Op2 sa 0) universe. 



Table 3: Cosmological parameters for various Op fitting \g\ < 0.1 



Op 


q 


P 




Oa 


HT 


w 


l5l 


0.3 


-0.539 ±0.011 


-0.705 


0.024 


0.696 


1.875 


-0.693 


0.052 


0.2 


-0.694 ±0.006 


-0.806 


0.020 


0.798 


2.796 


-0.796 


0.048 


0.1 


-0.847 ±0.003 


-0.909 


0.020 


0.899 


5.436 


-0.898 


0.054 



In the same way, for every given q there is a range of values for Op that satisfy the inequality (|62(l . As an example 
consider q = -0.3 and 7 = 0. From (jHlJ) we get 0.441 < Op < 0.5, while from ((23) we have 0.358 < Op < 0.466. 
Therefore, the physical range reduces to 0.441 < Op < 0.466. Consequently, for a fixed q the upper limit of Op is 
imposed by p > 0, while the lower limit by \g\ < 0.1. 

In table 4 we use p > and \g\ < 0.1 to constraint Op for different values of q. Once more, for a given q the values 
of Op are spread over a very small range, so the confidence interval for the density parameter is up to 97% — 98%. 

In order to put the above results in perspective, let us keep in mind that in FRW cosmologies q and Op are 
connected through the relation q — Op(l + 7)/2. Which means that Op (or q) can be determined by measurements 
of q (or Op) alone - in principle nothing else is needed. Unlike this, in our model q and Op are, from a mathematical 
point of view, independent from each other. However, from a physical point of view, they are not autonomous. The 
positive-energy condition creates a strong linkage between them, namely (|54|l . which is further strengthen by the 
observational evidence that \g\ < 0.1. What is amazing here is that these physical restrictions are strong enough 
as to allow the determination of one of these quantities by measurements of the other one, with a high degree of 
precision which today is close to 100%. 



Table 4: Cosmological parameters for various q fitting \g 


< 0.1 


q 


Op 


(3 




Oa 


HT 


w 


l5l 


-0.3 


0.454 ± 0.013 


-0.547 


0.028 


0.540 


1.252 


-0.533 


0.046 


-0.5 


0.326 ±0.007 


-0.678 


0.022 


0.670 


1.743 


-0.666 


0.046 


-0.7 


0.196 ±0.004 


-0.811 


0.020 


0.802 


2.847 


-0.800 


0.050 



*A comprehensive and updated discussion of the various experimental and observational constraints on the value of g (as well as on 
the variation of other fundamental "constants" of nature) has recently been provided by Uzan 1441 
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To finish this section we would hke to use the above results to "improve" the limits discussed in (|59|l and (jHOJ. 
From table 3 we see that the highest q, corresponding to Hp = 0.3, which satisfies the condition \g\ < 0.1 is q ^ —0.528, 
not^ —0.448 as in JHOJ. Since the observational data suggest that 0.3 < < 0.1, it follows that q = —0.528 is the 
upper limit for the acceleration parameter. Thus, 

g = -0.614 ±0.086. (64) 

Similarly for the density parameter. From table 4 we see that the lowest value of Hp, for q = —0.7, which satisfies 
(|62|l is rip = 0.192. The models for q = —0.7 with flp E (0.180, 0.192) are excluded because they do not comply with 
\g\ < 0.1. Consequently, from —0.7 < q < —0.3 we get 

0.192 <np< 0.3. (65) 

Consequently, our model together with the observational fact that |g| < 0.1, allows us to confine the whole range of 
possible values for the parameters q and Qp with and accuracy of about 86% and 78%, respectively. 

7 Past, present and future evolution 

The aim of this section is to outline the evolution of the universe according to our solution. 



7.1 The past 



In braneworld models, the early universe is dominated by the quadratic correction term. Accordingly, for t 
the "very" early universe, our solution (j25|l simplifies to 



ait) 



(7 + 1) Dt 



l/3(/3+l) 



Thus, 



q W 2 + 3/3, i? : 



1 



1 



The matter energy density and the vacuum energy density are given by 

2 _ 2(7-/3) 



Consequently, the cosmological term A(4) and gravitational coupling G vary as 

3(7-/3)' 



A 



(4) 



(7 + 1)' 



8ttG 



fc^5)(7-/3) 



7 + 1 



H. 



In addition. 



np = 



2(/3+l)(7-/3) 
(7 + 1)' ■ 



0, in 



(66) 



(67) 



(68) 



(69) 



(70) 



We note that in this limit, the solution is indistinguishable from the one for A(5) = (Tq = presented in (|23|l and 
(|24|l . This is physically reasonable because at this stage the evolution is dominated by p^, while the effects of the 
vacuum energy become important later. As an illustration, let us take a "small" Hp (say Hp = 0.01). Then, from 
(|7n)l it follows that /3 « 7 - Hp{l + -f)/2. Consequently, p/a = (/? + l)/(7 - /?) « 2/Hp « 200. 

^The rest of the values of g G (-0.528, -0.448), allowed by p > 0, in table 1 for flp 
with 15^ . 



0.3, are excluded because they do not comply 
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For f]p = 1/2, it follows from ^ that (3 
models, viz.. 



(7 — l)/2. Consequently, the scale factor matches the one in FRW 

a(i)^i2/3(7+l), (71) 

Hence, the deceleration parameter (|67ll reduces to q w (l+37)/2 which is the usual one in FRW cosmologies. However, 
in this limit neither the gravitational coupling is constant, nor the cosmological term becomes zero. Instead we have 



8ttG = -^H, 



A 



(4) 



4^ ' 



H 



1 



3(7+l)t' 



(72) 



where H is the Hubble parameter in FRW models. We note that this solution was also found in our recent work |34j . 
It describes the stage of the evolution when there is a perfect balance between ordinary matter and vacuum energy, 
namely p = a. At this stage there is no difference between the equation of state for ordinary matter and the one for 
the whole universe because w = 7. 

Throughout the early evolution, p > a, the parameters q, 7 and Qp satisfy (|56|l . The expansion {H > 0) is 
slowing down (q > 0) due to the huge gravitational attraction produced by the p^-term in the generalized Friedmann 
equation. As the universe expands; flp2 and (consequently) q decrease. At the same time ftp and Qj^ increase, such 
that flp2 + flp + = 1. However, flp is bounded above^, while flA is not. The upper limit flp = 1/2 corresponds 
to p = a. The epoch of (decelerated) expansion does not finish here; it continues until f2p « JIa, which denotes 
the moment when the universe becomes vacuum dominated^ In this phase of evolution the equation of state of the 
universe becomes quintessence-like {w < 0). 



Table 5: Outline of the evolution of the universe 


np 




q 


f3 






(i + 7)i/r 


w 


0.01 


2.985 ±0.000 


1.985 


-0.363 X 10-' 


198.00 


0.363 X 10--^ 


0.334 


0.990 


0.1 


2.846 ±0.004 


1.846 


-0.015 


17.973 


0.001 


0.345 


0.897 


0.3 


2.499 ±0.050 


1.499 


-0.063 


4.553 


0.017 


0.378 


0.666 


0.48 


2.04 ±0.24 


1.04 


-0.170 


1.833 


0.080 


0.446 


0.360 


0.490 


1.988 ±0.276 


0.988 


-0.191 


1.705 


0.092 


0.458 


0.325 


0.499 


1.909 ±0.342 


0.909 


-0.230 


1.550 


0.114 


0.479 


0.273 


0.5 


1.5 ±0.75 


0.5 


-0.5 


1 


0.25 


0.667 





0.499 


1.091 ± 0.342 


0.091 


-0.502 


0.458 


0.387 


0.797 


-0.273 


0.490 


1.013 ±0.274 


0.013 


-0.516 


0.378 


0.417 


0.846 


-0.325 


0.48 


0.960 ±0.240 


-0.040 


-0.530 


0.333 


0.440 


0.886 


-0.360 


0.3 


0.461 ± 0.011 


-0.539 


-0.705 


0.024 


0.696 


1.875 


-0.693 



The quantities p/a and J^a do not depend on the specific value of 7. Unlike this, q, (3, w and T do depend on it. In 
table 5 we present their values calculated for 7 = 0. For any other 7 the corresponding parameters q-y, (3-y and 
are given by 

g^ = g(l + 7) + 7, = /?(1 + 7) + 7, w^=w;(l + 7)+7. (73) 



7.2 The present 

Current dynamical mass measurements suggest that the matter content of the universe adds up to 30% of the critical 
density^. Correspondingly, if we take Vip = 0.3, then the universe is speeding up its expansion with q = — 0.539±0.011 
(see table 3), which is consistent with recent observations. 

Although the accelerated expansion of the universe, at the current epoch, is driven by the dominating dark energy 
associated with A(4), the dynamical evolution is established by the total matter content, which includes regular (cold 

®The models for which 1/2 < Qp < 1 never become vacuum dominated. These are of no interest for us here. 

^We note that cr > p, does not necessarily imply f^A > f^p- This is a consequence of the quadratic correction term in the Friedmann 
equation. 

^Radiation 0.005%, ordinary luminous baryonic matter 0.5%, ordinary non-luminous baryonic matter 3.5% and exotic (non-baryonic) 
dark matter "observed" through their gravitational effects 26%. 
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dark plus baryonic) matter as well as the dark energy component. The effective equation of state of the universe 
for q = -0.539 ± 0.011 is restricted to be 



- 0.7 < i^yj < -0.685. (74) 

Going aside for a moment, let us to note that in our model the equation of state of the universe, pt = wpt, is similar to 
the one derived for quintessence in the framework of flat FRW cosmologies j^Hj- In this framework, observations from 
SNe la and CMB indicate that the equation of state for quintessence wq = VqI PQ has an upper limit wq k. —0.6, 
which is close to the lower limit wq « —0.7 allowed for quintessence tracker fields [231 • Since the dynamics of 
spacetime is governed by the effective equation of state — wpt (not only by the vacuum energy), the fact that 
Wq « u; is a clear evidence of how little is the effect of ordinary matter on the present evolution of the universe. 

For practical reasons, it is convenient to express the solution in terms of the quantity z, the redshift of the light 
emitted at time t, which is the directly observed quantity used in research on distant objects. The connection between 
the scale and the redshift parameter is given by 

«(^) = (75) 



where t is the present time, i.e., t — T, so that z — now. From H25|l and (|30|l we get 

^'+^^=7^^h^' T^fc2 ao(/3 + l)t, (76) 



where for convenience we have introduced the constant i] as 

477 = f 2 + 4f. (77) 
Employing the assumption 7 = today, and (|75(l , the cosmological parameters in terms of z are 



6 



2 V [?7+(l + z)3(i+«]y ' 

l + t^(L) (1 + ;,)3(l+/3) 



13 + 1 \ao 



87rG(z) 

^""^'^ ~ [2(/3+l) + (l + z)3(i+/3)] 

n^iz) _ {fj + l)G{z) 



/3[2(/3+l)-/3(p/ao)(l + z)3(i+/3)] 

(78) 



np [f^+(l + ;2)3(l+/3)]G' 

where p , G and flp represent the present values oi p , G and Op, respectively. Thus, in order to be able to predict 
the above parameters, at any z, in terms of their present values, we need to know the constants A:(5), [3 and ctq- 

7.2.1 Evaluation of the constants 

In order to evaluate the various constants as well as the vacuum energy and the cosmological "constant" today, let 
us take Clp, q and (3 from the last row in table 5. For these values, our universe is almost three times older than the 
dust FRW universe, for which HTpRw — 2/3. Besides, the vacuum energy density is much bigger that the matter 
energy density, namely, 

CT«41.67p, p = 3.581 X 10"^ (^) . (79) 



. G 

The five-dimensional quantities fc(5) and A(-5-) can be evaluated as 



kf^) = 0.363 = 10.146 , A(5) = -5.390 x lO^{Hp) = -19.303 ■ (80) 
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If we take G = 6.67 x IQ-^ crn^ gr'^ s''^ and iJ w 0.7 x IQ-^^yr-^ ^ 2.22 x 10~i*^s-\ then 

p = 2.65 X IQ-^^grcm-^ 1.14 x IQ-^'^GeV^. 
The universe is around 26.8Byr old and, and the rest of the parameters are as foUows 



(81) 



A(4) 


= 1.03 X lO^^-'^s"^ = 


4.46 X lO^^^Gey^ 


PA 


^ 6.14 X lO-^°grcm- 


-3 2.65 X IQ-^'^GeV^, 


a 


= 1.08 X IQ-^^grcm' 


= 4.66 X IQ-'^^GeV'^, 




^ 1.02 X IQ-'^^grcm^ 


-3 = 4.40 X IQ-^^GeV^, 


'^(5) 


= 3.05 X 10"cTO^gr" 


^s-^ = 4.66 X lO^Gey-^ 


A(5) 


= -3.17 X lO^^'^cm" 


^grs-^ = -0.90 x lO^^^GeF^ 




= -9.67 X lO^^'^s"^ 


= -4.19 X IQ-^^GeV^. 



(82) 

In table 6 we illustrate how a modest change of q, within a range of allowed values, influences some of the parameters 
in the solution. 



Table 6: Dust universe for ftp = 0.3, and various q 



q 


/3 






HT 


w 


l5l 


-0.536 


-0.707 


0.031 


0.695 


1.843 


-0.691 


0.066 


-0.543 


-0.703 


0.015 


0.698 


1.932 


-0.695 


0.033 


-0.550 


-0.700 


0.000 


0.700 


2.222 


-0.700 


0.000 



Turning back to the formulae in terms of the redshift we can use ()82(l to obtain the constants in H78|l . In 
particular, for fj we obtain 



22.7715, 



(83) 



and {p/o-q) — 0.02596. With these values, we can now use H78II to reconstruct the past of the universe. This model 
has a number of interesting properties but we leave their discussion to another opportunity. 

7.3 The future 

For t ^ to, far from the initial singularity, a ~ ao and a ^ p. The asymptotic evolution of the expansion scale 
factor, as i — > oo, is given by 

a(i)-consixi2/3(/3+i)^ (84) 

irrespective of the equation of state of the matter on the branc. Thus, for the future evolution of the universe, we 
obtain a cosmological model with constant G and a non-vanishing cosmological term A(4), which varies with time. 
Namely, 



87rG 



r,4 f^O 



6 



A 



(4) 



3(7-/3) 
(7 + 1) 



H\ H 



1 



3(/3+ 1) f 



The matter density and pressure become 

8nGp 



3(7 + l)(/9+l)t2' 



P = IP- 



(85) 



(86) 



Here /3 is related to the deceleration parameter as /3 w (2g — l)/3. We also note that the ratio G/G varies from 
(G/G) w -3(/3 + 1) for t « to (G/G) ^ for < ^ oo. 

At this stage, the parameter /3 becomes totally determined by f^A- Namely, for any given 7 



(5 = — ilA = w. 



(87) 
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The corresponding effective, or total, energy density is obtained from (|85|l and H8t)|) as follows 



%-KGpt = A(4) + 87rGp ^ 3(^,^1)2^2 ' = '^^Pt- (88) 

In this limit fip + ^A = 1, which means that the brane-model under consideration grows up to become, asymptotically 
in time, indistinguishable from a spatially flat FRW cosmological model based on a mixture of cold dark matter and a 
variable cosmological "constant" . The model can be completely determined by measurements of only one cosmological 
parameter. For example, if we know Jlp, then from 1)51(1 and ((52|l we get q and w, respectively. This in turn gives /3 
from ((87() . which allows to calculate the age of the universe and the Hubble parameter. 



8 Summary and conclusions 

An important feature of braneworld theories is that the gravitational coupling G and the effective cosmological term 
A(4) are related to tr, the vacuum energy density of the 3-brane. They are not independent, as in Jordan-Brans-Dicke 
and other multidimensional theories Thus, in brane-theory either all of them are constants or they vary 

simultaneously. 

In the case where the vacuum energy is constant, and the ordinary matter satisfies the equation of state p = 7p, 
the conservation equation Tj^.^ = leads to p a^'^('''+^'. In this paper, in order to integrate the generalized FRW 
equation of braneworld models, we have assumed p ^ a^^^P^^\ where /3 ^ 7 is a constant. We have seen that /3 « 7 
at the early universe, while /3 « — f^A at late times. This assumption constitutes the "minimum" extension to the 
usual braneworld models with constant vacuum energy (/3 = 7). It adds new algebraic possibilities, while keeping 
the calculations simple enough as to lead to physical effects that can be studied exactly. 

The expansion of the universe at late times is determined by the constant Ag/ / introduced in (|21|l . We have set 
Ag/j = 0, thus excluding from the outset the existence of recoUapsing or exponentially expanding solutions, and 
studied in detail the spatially-flat cosmological model with C = 0. 

In general is a "function" of the cosmological parameters q and Vlp. If /5 = 7, then from H4()|l we obtain 
q = 2 — 3flp/2, in the dust dominated era. Thus fixing Qp today also fixes q. Notice that q is positive for any physical 
value of Sip, meaning that a brane-universe with constant vacuum energy must be slowing down its expansion. 
However, for /3 7^ 7, this is no longer so; the cosmological parameter J^a is now a dynamical quantity which changes 
this picture. The acceleration parameter is now given by H51|) and (|73|l . so that as the universe expands q decreases 
from q-y = 87 + 2 at the beginning of the universe to q —^ —1 -\- 3Qp/2 at late times. Thus, leading to a phase where 
the universe, instead of slowing down, is speeding up its expansion. 

As a consequence of the introduction of /3 7^ 7, the cosmological parameters q and Qp are independent from each 
other. However, for any given value of flp (or q), the positivity of energy and the observational limits on the variation 
of G, allow to determine q (or Hp) with a precision which is close to 100%. 

Therefore, the equation of state of the universe H52|l varies very slowly in each epoch of the universe (see table 
6). But it does change when 7 changes, from ((49(1 it follows that w is affected by 7 through its influence on q. 
As the universe expands the equation of state changes from w w 1 + 27 at the very early universe, to become 
quintessence-like at late times, namely, w ^ —1 + ftp. Recent data from type la supernovae in distant galaxies and 
the cosmic microwave background constrain the equation of state of the universe (35| . If the universe is assumed to 
be spatially flat, as we do here, then the upper limit for w is —0.6, which is close to the lower limit w « —0.7 allowed 
for quintessence tracker fields jSE]- These limits^ are quite similar to those derived here (|7^ . 

A remarkable feature of our model is that f2p is bounded above. Consequently, at some stage of the expansion 
the vacuum energy starts dominating the evolution. Precisely, the density parameter reaches its maximum value 
when the epoch where p > a comes to an end. From then on cr > p. One might think that at this stage something 
radical happens. In particular, that the evolution becomes vacuum dominated (I^a > ^p)- However, because of 
the quadratic correction which is still considerable at this stage, this is not so. The reality is that for most of the 
evolution of the universe f2p > J^a- Only recently, for HT « 0.8 (see table 5), the black energy overtakes the matter 
energy and drives the universe into a period of accelerated expansion. 

^ Since > 0, our braneworld model excludes phantom energy, which requires ui < — 1 
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The model discussed here grows to be FRW-like in time, in the sense that the quadratic correction Qp2 becomes 
negligible and G becomes constant at late times. However, there is an important distinction; the cosmological 
"constant" varies with time and is proportional to H^. It should be noted that the relation A(4) « const x is 
favored by observations and there is an extensive literature suggesting that it plays a fundamental role in cosmology 

mm 

Another important feature of the model is that it contains no "adjustable" parameters. All the quantities, even 
the five-dimensional ones, can be evaluated by means of measurements in AD. We illustrate the values for fc(5) and 
A(5) in l|8()|l and H82|l . We remark that these are "universal" constants fixed by the five-dimensional embedding bulk. 
Regarding the "cosmic coincidence" problem (TO!, it is interesting to note from H82|l . that in our model the missing 
energy density pA is close to p, the matter density today, without any (obvious) fine-tuning of their ratio at the early 
universe. 

Let us briefly compare this work with our previous investigation 34 . In that paper our working hypothesis was 
that [G/HG) = g, remained constant throughout the evolution of the universe. Then, the observational constraint 
\g\ < 0.1 required ctq = 0. This means that, as the universe expands a —> and, from H12|) . A^^^ k'^^^A(^^/2. 
If our universe is assumed to be embedded in an anti-de Sitter five-dimensional bulk, then A(-4') changes its sign, 
from positive at the beginning to negative later, which in turn causes the universe to recoUapse in the remote future 
(HT « 10.76). For a universe embedded in a bulk with A(5) = 0, we obtained the ever-expanding solution given by 
(j23|l and (|24|l . In the present work ctq 7^ and there is neither recoUapse, nor exponential expansion. We excluded 
both by setting A^jf — 0. 

In spite of this, and the fact that the recoUapsing model is few times older than the present one, the cosmological 
parameters q, J^a, f2p2 as well as the five-dimensional quantities fc(5) and Ajs) have similar values in both models. 
Therefore, we cannot discriminate between them on the basis of the cosmological parameters only. In order to do 
that, another type of analysis should be done, for example the behavior of perturbations and structure formation. 
These problems are beyond the scope of the present work. 

In summary, in this paper we have constructed a cosmological model, in the context of five-dimensional braneworld 
theory, where the vacuum energy varies with time. The model provides a simple scenario where we can study the 
simultaneous variation of G and A(4) in an universe which is spatially flat and expanding with acceleration, in 
agreement with recent observations. The vacuum density a may vary in models where (744 = — <I>^ is a function of 
time. In this connection, it is important to mention that the ratio {^/^) appears in different contexts, notably in 
expressions concerning the variation of rest mass and electric charge. Regarding the time- variation of the gravitational 
"constant" G, in different models with extra dimensions the ratio [G/G) is found to be proportional to (^/^). 

Therefore, we have a scenario where the observed acceleration of the universe is just one piece in the dynamical 
evolution of an universe where the so-called fundamental "constants" are evolving in time. The fact that asymptoti- 
cally in time G becomes constant, suggests that the size of the extra dimension stabilizes in time and, consequently, 
the rest mass and the electric charge become constant too. This study can be extended by considering other pos- 
sible functional forms for the vacuum energy, obtained from realistic potentials V{(j)) determined from the observed 
cosmological functions. 



Appendix: The vacuum energy as the energy of a slowly-evolving scalar 
field 

Although the energy-momentum tensor for vacuum looks like an ideal ffuid with negative pressure, the dynamics 
of a varying vacuum cannot be that of a fluid because stability would require dp/dp > 0. The simplest working 
model for a dynamical vacuum energy is the dark energy of a single homogeneous scalar field (j) (its slight clumping 
is neglected) with self-interaction potential V{(j)). 

For completeness, we would like to show here that the vacuum energy density H17(l can be derived from a simple 
potential. 

For the cosmological metric Q the equation governing (j) is 

0-^3^0+^ = 0. (A-1) 
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The energy density and pressure are given by 



1 



0' - V{^). 



(A-2) 



Here the appropriate dimensions for (j) are cmT^/'^gr^/'^s^. If V^(0) is sufficiently flat, then will evolve very slowly 
and (jP' ^ V . In this scenario the field energy approximates the effect of a cosmological term with « 
Combining the above equations we get 

J.2 



Consequently, 



V = 



a dp^ 
3 da 

^ a 
6 da 



We now identify the vacuum energy with the energy of the scalar field, viz., a 

(7-/?)(l-/3) 
V ^ aQ + ^ p. 



2(/3+l) 



(7 - P)P- 



For the solution under consideration (|25|l we get 



(7-/3) (1-/3)^0 
(7+l)sinh2C(0-0o)' 



(A-3) 
(A-4) 

Then using p7|) we obtain 

(A-5) 
(A-6) 



where 0o is a constant of integration and^'' 



fc45)(To(/3+l)(7 + l) 



/ 8(7 - P) 

The time evolution of the scalar field is given by 



0.704 X IQ-^c-nr'/'^gr 



(r + 2) + Vt2 + 4r 



(A-7) 



(A- 



for most of the 



We notice that V{4)) sharply decreases with the increase of 0, so that V{4)) « o-q and « 
evolution of the universe. Potential ljA-6p takes a more familiar form in the early universe and the late universe. 

Power law expansion In an early universe, dominated by matter or radiation, the term drives power law 
expansion a ^ 

ii/3(i+/3). At late times the dark energy becomes dominant and also drives a power law expansion, 
viz., a ^ t2/3(i+/3)^ general, for a scale factor a ^ t"" the field equation l|A-ip has two types of solutions. The first 
one corresponds to the pure inverse power-law potential, namely, 



const 



t 



2/(2+q) 



(A-9) 



where a is a dimensionless parameter related to the equation of state of the universe. The second solution is for the 
exponential potential 

m ^ .0 + fc/^l^^^e-^^ C(^-0o)=lnr. (A-10) 
(7 + 1) 



-'For the evaluation of C we have used the numerical values found in 1821 . 
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In our solution, the potential ljA-6|l can be expanded in powers of (f>. The first term of the expansion corresponds to 
(|A-9p with a = 2. On the other hand (|A-10p is the asymptotic (r — > oo) form of (|A-6p . 

Therefore, for the largest part of the evolution the scalar potential (|A-6|) is a good model for the cosniological 
term. But, notice that it does not work at the very early universe where « —/3p<p ~ ~1P4>'-! another mechanism 
for a variable cosmological term should be used for that era. 

Before finishing, we would like to mention an important difference^^ between the present model and the FRW 
cosmologies with a variable cosmological term, modeled as a scalar field. Namely, in FRW models A(4) is a linear 
function of the vacuum energy, viz., A(4) ~ p^, while here it is a quadratic function; A(4) = k^^^{a'^ — CTo)/12, with 
a = p,p. In spite of this, the numerical values for the vacuum energy and cosmological constant as calculated in H82|l 
are similar to those expected from other branches of physics. 
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